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1. Introduction 

Although there is strong evidence that Sr2Ru04 is a triplet superconductor 
[2, 3], the full symmetry of the equilibrium state below Tc remains open 
to debate [2-14]. There exist strong indications for broken time-reversal 
symmetry in the superconducting state [15, 10] and equally convincing 
measurements showing that the order parameter d(k) has a line of nodes 
on the Fermi surface [16-19]. The reason why this state of affairs represents 
a puzzle is that for all odd parity spin triplet pairing states in tetragonal 
crystals, group theory does not require the simultaneous presence of both 
broken time-reversal symmetry and line nodes [20]. To explain this incon- 
sistency different three dimensional models of pairing have been proposed 
[11-13]. In fact the experimental results on heat transport [19] seem to 
favour the horizontal, with respect to (ab) crystal plane, line nodes. Usu- 
ally in a multi-band BCS like model, with different coupling constants for 
each band, one generically finds multiple phase transitions as the different 
sheets of the Fermi surface are gaped on lowering the temperature. Since 
experimentally there is only one jump in the specific heat, at Tc = 1.5 K, in 
constructing a sensible model one must eliminate such multiple transitions. 
Zhitomirsky and Rice [12], in their simplified two band model, considered 
an interaction which couple the order parameters of the different symmetry. 
The presence of this inter-band interaction leads to a single superconduct- 
ing transition. Microscopically such coupling of different bands comes from 
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effective three site interactions. They assumed that one band in Sr2Ru04 
is the most important for superconducting pairing, while two other are 
gaped via an inter-band proximity effect. On the other hand Annett et 
al. [13] have proposed sightly different model. They considered a three 
orbital three dimensional model with effective in plane and out of plane 
nearest neighbour interactions, but did not allow for interactions mixing 
the symmetry of the order parameters. The price to pay in this model is 
the fine-tuning of two interactions in order to have a single superconducting 
transition at Tc ~ 1.5 K. The calculations show that the model [13] explains 
quantitatively the T dependence of the specific heat, penetration depth and 
thermal conductivity without additional fitting parameters. 

It is the purpose of this work to extend the model of Annett et al. [13] 
by allowing for symmetry mixing interactions. To this end we add small 
three point interaction to describe possible inter-orbital proximity effects. 




Figure 1. Calculated tight binding Fermi surface of Sr2Ru04. 



2. The orbital model of superconductivity 

To describe the superconducting state we use a simple multi-band attractive 
Hubbard model. Its Hamiltonian consists of two parts: 

H = Ho + Hint, (1) 

where Hq denotes the tight binding electron part corresponding to the 
experimentally observed band structure of Sr2Ru04. 

^^0 = ii^m- fJ')Sijdmm' -tmm'{ij))Ci'„ia^jm'a, (2) 

ijrnrn' .(J 

here m, m' = a,b,c refer to the three ruthenium t2g orbitals a = xz, h = yz 
and c = xy and i and j label the sites of a body centered tetragonal lattice. 
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Hma ^i'm.a ^re the Fermion creation and annihilation operators for an 
electron on site i and orbital m with spin a. 

The hopping integrals tmm'iij) and site energies Em were fitted to re- 
produce the experimentally determined Fermi surface [21]. The calculated 
Fermi surface is shown in Fig. 1. 
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Figure 2. Body centered tetragonal lattice of Ru atoms in Sr2Ru04. Sites i, j, I 
correspond of possible realization of three nearest lattice sites. Note that two of sites 
must lie in Ru02 plane, while the third lies outside it. 

The second part H-int (Eq. 1) describes the general electron-electron 
interaction and can be written as 

Hint = — 2 X] ^nun'oo'{'''j^^)'4ma^jrn'actoa'^so'a'- (3) 
ijlsmm'oo'cra' 

Some of the interaction constants U^^,^^,{ijls) are assumed to be at- 
tractive: those acting between electrons on nearest neighbour sites with 
parallel spins and responsible for superconductivity. In Eq. 3 all site indices 
ijls can correspond, in general, to different sites but we assume, as is usually 
the case in an isotropic substance, that n-point interactions satisfy the 
relation 

|C/(n=l)| > \U^ri=2)^ ^ |C/("=3)| > \U^'^=^)\. (4) 

A single point interaction {U^^^) refers, as usual, to the on-site repulsive 
Coulomb repulsion which does not contribute to electron pairing while the 
U^"^^ , assumed to be attractive, represent the main contribution. The three 
point interaction coupling (U^'^^) is apparently smaller than the bond term 
(U^'^^), however its influence on pairing mechanism may be important [12], 
so we retain it, neglecting a four point coupling U^^\ 

Thus, the interaction part of Hamiltonian in our case reads: 

Hint ~ 2 ^ ^ Ujy^^iiij^flima'f^jm'a' 
ijmm'acr' 

ijlmm'aa' 
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where U!^^,{ij) (for i ^ j) describe an attraction between electrons on the 
nearest sites with spins a and a' and in orbitals m and m' while U!^^, (ijl) 
{i 3 constant is the interaction (of any sign) between three nearest 
neighbour sites respectively (Fig. 2). It is called the assisted hopping term 
[22]. In our model, as will be disscused later, we assume that this term 
couples the in- and out-of- plane order parameters. 

The actual calculations consist of solving, self-consistently, the following 
Bogoliubov-de Gennes equation: 



(TIT 



where H^^i{ij) is the normal spin independent part of the Hamiltonian, 
and the A^^, (ij) is self consistently given in terms of the pairing amplitude, 
or order parameter, A^^,(ij) is defined by the usual relation: 

+ J2Umm',oo'mXoo'iil), (7) 



where 



and 



loo' 



x^'^'m = Y.<movrm'A^ - vm) , (§) 



is Fermi function, /3 = l/ksT, Hb is Boltzmann constant and v enumerates 
the solutions of Eq. 6. 



3. Symmetry of order parameters 



We solved the above system of Bogoliubov de Gennes equations (6-8) 
including all three bands and the experimental three dimensional Fermi 
surface (Fig. 1). We assumed that the pairing interaction U'^^,{ij) for 
nearest neighbours in plane is only acting for the c {dxy) Ru orbitals and 
the nearest neighbour inter plane interaction acts only in a and h orbitals 
{dxz, dyz)- The motivation for this is that the dominant hopping integrals in 
plane are between c orbitals, and the largest out of plane hopping integrals 
are for a and h. On the other hand the additional three point, assisted 
hopping, interaction provides coupling between different orbitals a,b and c 
Therefore we have only three coupling constants U\\ and Ux_ describ- 
ing these physically different interactions in- and off-plane and Ui which 
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correspond to inter-orbital coupling. Our strategy is to adjust these phe- 
nomenological parameters in order to obtain one transition at the exper- 
imentally determined Tc- Of course, this can be done for many choices of 
interactions parameters. Compare the results of our previous orbital model 
[13] with these obtained here for generalized orbital model with three point 
interactions taken into account. 
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Figure 3. Order parameters, |A^„| = |A^i,| and |AJc| as functions of temperature for 
various sets of interactions: U±, U^, Ui ((a) for 0.590t, 0.494f, 0.0; (b) for 0.400, 0.493t, 
0.004t; (c) for 0.400t, 0.400t, 0.054t; respectively; where t = 0.0816 eV). 

Because the pairing interactions U^^^'^^,{ij) were assumed to act only for 
nearest neighbour sites in or out of plane, the pairing potential A^^,(ij) 
is also restricted to nearest neighbours. We further focus on only odd 

parity (spin triplet) pairing states for which the vector d (0,0, d^), i.e. 
/(il) = A''^^ Aij), and A^^ = A^^ ,(ij) = 0. Therefore in general 

mm V / mm V / ' mm \ •J / mm \ •J / o 

we have the following fourteen non-zero order parameters (i) for in plane 
bonds: Acc{ex), Acc{ey), and (ii) for inter-plane bonds: AaaCRij), AabiRij), 
AbbCRij) for Rij = (±a/2,±a/2,±c/2). 

Taking the lattice Fourier transform of Eq. 7 the corresponding p-wave 
pairing potentials in k-space have the general form of order parameter 
(suppressing the spin indices for clarity): 

Acc(k) = A^^ sin + A^^ sin ky (10) 
for the c orbitals and, 

Amm' (k) = (^A^^, sin y cos ^ + A^^, sin ^ cos cos ^ (11) 

for m, m' = a,b. A more complete k-space analysis can be found in Ap- 
pendix A. 

Note that beyond the usual p-wave symmetry of the sinkx and sinA;^ 
type for the c orbitals, we include all three additional p-wavc symmetries 
of the sin/cj/2 type which are induced by the effective attractive interac- 
tions between carriers on the neighboring out-of-plane Ru orbitals. These 
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Figure 4- Eigenvalues on the a (a), l3 (b) and 7 (c) Fermi surface sheets, E(kF) along 

kz direction for interactions chosen as in Fig. 3a, (d) ECkp) for /3 band Fermi surface 
and interactions chosen as in Fig. 3c (T = 0) . 



interactions arc also responsible for the f- and pz-wave symmetry order 
parameters [13]. These relatively small components with smaller Tc will be 
neglected in this paper. The motivation comes from the strong impurity 
effects in the material [23]. 

We have solved Eqs. 6-8 for various sets of interactions: U±, ?7||, Uj. For 
Uj = the order parameters have the symmetries A^^ = iA^^, Aff^ = iA^^ 
as expected for a pairing symmetry [4] of (A;,^. + iky)ez type corresponding 
to the same time reversal broken pairing state as '^He-A. The off-diagonal 
components, such as A^^ are small but non-zero, as are A^^ and A^„. Note 
that the k-space pairing potentials Ajnm'O^) do not directly correspond 
to the energy gaps on the Fermi surface sheets, because the tight-binding 
Hamiltonian is non-diagonal in the orbital indices. Instead in Fig 4a-c we 
show the energy eigenvalues of our Hamiltonian evaluated at the Fermi 
surface for a, [3 and 7 bands. One can see that there are horizontal circles 
around the cylindrical [3 Fermi surface sheet (Fig. 4b), while the order 
parameter at the 7 sheet (Fig. 4c) is node-less. Interestingly, in our model 
a sheet does not posses line nodes. Due to the small diagonal distortion of 
a Fermi surface (Fig. 1) it has point nodes (Fig. 4a). 

^11 ) Ui were chosen in such way as to give a single transition with 
Tc = 1.5 K. Figure 3a shows various A obtained for our [13] orbital model 
(JJj_ = 0.590t, C/|| = 0.494t but Ui = 0.0). Here the interactions were tuned 
separately to give approximately the same Tc- All components of the order 
parameter in a,b and c orbitals are active. 
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Figure 5. Temperature dependence of (f> (A^c — e"^Acc). Interaction parameters U±, 
U\\, Ui as in Fig. 3 for (a), (b) and (c); respectively. Note that for Ui = (a) 4> = -k/2 
and AK, = iA^,. 



In the case of Fig. 3b U±, 17^, Ui are equal to 0.400f, 0.493t, 0.004f, 
respectively. Here a,b orbital order parameters A^^, A^^ have much smaller 
values than A^^ and only due to small three point inter orbital coupling 
Uj. We have the proximity effect and hence a unique Tc- 

In the last case C/^=0.400f, U\\ = 0.400t, Uj = 0.054t (Fig. 3c) the 
situation is different. To get Tc = l.bK both interactions U±, C/y take on 
smaller values comparing to those in Fig. 3a. The increase of the inter- 
orbital interaction [// changes Tc and the values of A^^, (T = 0) as well 
as their temperature dependence. The order parameters A^„ fulfill more 
general relations A^^ = e"'^A^^ and A^j^ = e"^AS(^, where 4>\s a. temperature 
dependent phase. This dependence for all three cases are plotted in Fig. 5. 
For [7/ = (/> is constant and equal to 7r/2 (curve a). A similar behaviour 
with « 7r/2 can be found for small Ui (curve b) while for relatively 
large Ui {Uj = 0.054t) depends strongly on temperature (curve c). It 
is worthwhile to note that for low temperature T all order parameters are 
complex. It is also true for higher T but for Ui equal to or relatively small. 
For larger interaction [// there is a temperature T* (T* < Tc) above which 
the order parameters are real. This is a source of additional vertical line 
nodes appearing at the Fermi surface. For Uj = 0.054t the horizontal line 
nodes are still present on the cylindrical /? Fermi surface sheet (Fig. 4d). 
Due to (j) / 7r/2 the eigenvalues on the Fermi surface show a two-fold 
symmetry instead of four- fold one present for Uj = (Fig. 4b). 

To decide which of these three cases is closer to experiment we calculated 
the specific heat for the same interaction parameters as above using the 
relation: 

C=-2kBP'Y.Em,.^^^§^. (12) 

In Figs. 6 a-c we have plotted the specific heat versus temperature. In 
all cases the low temperature limit of the specific heat is power law, because 
our gap parameters have line nodes on /3 Fermi surface sheet. 
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Figure 6. Lines: Specific heat calculated for various sets of interactions (parameters 
U^_, U\\, Ui as in Fig. 3); full circles denote experimental data [16] 



Going to larger temperatures we observe in Fig. 6a, that C has a dis- 
continuity at Tc of 27mJ/mol, which is in agreement with the experiment 
[16, 13]. The slope of C/T versus T also comes out correctly. The values 
of the interaction parameters U± and U\\ arc in this model tuned to give 
single superconducting transition at temperature equal to 1.5 K. Interest- 
ingly, the two other figures (Fig. 6b-c) which correspond to the inter-orbital 
proximities have too small value of the specific heat jump at Tc- Technically 
it is due to different curvatures of temperatures dependences of ^mm'- 



4. Summary and conclusions 



In summary we would like to emphasize two points. Firstly, wc have ana- 
lyzed the orbital model of superconductivity with and without 'inter-orbital 
proximity effect'. Both versions of the model allow for complex order param- 
eter with line nodes on some of the Fermi surface sheets. Our description is 
a real space one with two point tight-binding interaction such as naturally 
arises in multi-band, extended, negative U Hubbard model [24]. 

Secondly, we wish to stress that in our original approach to the problem 
the parameters which describe the normal state are determined by fitting 
them to the very accurately known Fermi surface. The measured deter- 
mines both coupling constants J/y, U_\_. As for the physical mechanism of 
pairing, the fact that U±_ implies that the pairing interaction is fairly 
isotropic in spite of the layered structure of the system. The model with 
"inter-orbital proximity" effect requires even more fine tuning of parameters 
in order to get correct slope and jump of the specific heat. 

In conclusion, the scenario with active cc orbitals and "orbital proxim- 
ity" mechanism for superconductivity in a and b orbitals is not consistent 
with specific heat data on Sr2Ru04. 
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6. Appendix A 

In case of a pure system one can take a Fourier transform of Eqs. (6-8). 
Thus 

sr- ( E^- Hmm'O^) A-;,(k) \ / <^,,, \ 

^, 1, A-;*, (k) + (k) J [ ; ' ^ ■ ' 

The gap equation (Eqs. 7-8) can be rewritten in k-space as 

^mU^) = ^Ef^rm'(k-q)x;;l'(q) 

q 

+ ^ E Urnm',oo' (Q, k - ci)xZ' (q) • (A.2) 



N 

q,oo' 



and 



Xw(k) = U^maVt^'A^ - 2f{E^)) . (A.3) 



In the orthogonal crystal (Fig. 2) various matrix elements of the inter- 
action U responsible for p-wave paring can be written as 



C/,,(k,q) = 2[/||y(k)y(q) 

t^mm' (k,q) = 8[/^y(k)y(q) for m,m' = a,h (A.4) 

[/wcc(k,k') = 8[/jy(k)y(q) for m',m = a,6 

C4cmm'(k,q) = 8[7/y(k)y(q) for m',m = a,b 

where V(k.) and t^(k) can be expressed as 

y(k) = (sinfca; -I- sinky) (A.5) 

T>/i \ f ■ kx ky , ■ ky ky\ 

V[k) = ( sm -— cos -I- sm cos I cos -— . 

Thus, the general form of order parameter 

Acc(k) = sin + A^^ sin fej, (A.6) 
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for c orbitals and, 

(k) = (a^„, sin ^ cos I + A^^, sin ^ cos ^) cos ^ (A.7) 
for m, m' = a,b orbitals. 
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